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PART  A 

COMPUTATIONAL  STUDIES 
H.  M.  STERNBERG 


I .  INTRODUCTION 

Computational  work  with  Amatex-20  (20/40/40  RDX/TNT/AN)  was  under- 
taken in  connection  with  its  potential  use  as  an  explosive  fill  for 
the  BLU-61B  and  BLU-63B1 2  bomblets.  The  general  plan  of  the  work  was 
to  use  existing  cylinder  test  data  obtained  at  Picatinny  Arsenal7 
(PA)  and  the  Los  Alamos  Scientific  Laboratory3 * * *  (LASL)  to  get  an 
equation  of  state  scheme  to  use  with  Amatex-20  in  explosion  hydro- 
dynamic  calculations.  Computations  would  then  be  made  of  the  casing 
velocities  of  Amatex  20  loaded  BLU-61B  and  BLU-63B  bomblets. 

The  use  of  ammonium  nitrate  (AN)  raises  the  possibility  of  a so 
called  "non-ideal"  detonation,  i.e.,  part  of  the  energy  may  be 
released  after  the  detonation  front  has  passed.  In  many  applications 
this  late  energy  release  will  not  affect  the  conditions  at  the  deton- 
ation front.  It  can  be  shown,  however,  that  if  the  order  of  the 
reaction  is  equal  to  or  greater  than  one,  as  is  invariably  the  case, 
a steady  state  detonation  is  not  possible1*,  and  the  late  energy 
release  must  ultimately  affect  the  front  if  the  charge  diameter  is 
considered  to  be  infinite  and  the  detonation  is  allowed  to  propagate 
indefinitely . 


1.  Bomb,  Frag-BLU-61B,  Dwg.  No.  68D5900,  and  Bomb,  Fragmentation 
BLU-63B  or  BLU-86B,  Dwg.  No.  68D5994,  U.  S.  Air  Force,  Eglin  AFB, 
Florida . 

2.  J.  Hershkowitz  and  J.  Rigdon,  "Evaluation  by  a Modified  Cylinder 
Test  of  Metal  Acceleration  by  Non-Ideal  Explosives  containing  Ammonium 
Nitrate"  Picatinny  Arsenal,  Dover,  New  Jersey,  Technical  Report  4611, 
April  1974. 

3.  A.  W.  Campbell,  R.  P.  Engelke,  A.  Popolato,  and  T.  Rivera, 

"Physical  and  Processing  Characteristics  of  Non-Ideal  Explosives- 

Special  Emphasis  on  Amatex  (Task  A),  "Los  Alamos  Scientific  Laboratory, 
Los  Alamos,  New  Mexico,  LA-57  11-PR  Joint  Services  Explosives 

Program,  Quarterly  Report,  March  16  through  June  30,  1974,  compiled 

by  A.  Popolato. 
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Definitions  of  the  terms  "at  the  detonation  front"  and  "behind  the 
detonation  front"  are  needed  to  fix  the  ideas.  At  the  detonation 
front  will  refer  to  points  in  the  flow  where  the  particle  velocity 
plus  the  sound  speed  (u+c)  is  equal  to  or  greater  than  the  detonation 
velocity  D.  Behind  the  detonation  front  will  imply  u+c<D.  It  should 
be  borne  in  mind  that  the  distance  between  the  point  where  D=u+c  and 
the  detonation  front  will  vary  with  the  run  distance  (see  Ref.  4). 

For  practical  purposes  it  suffices  to  define  "at  the  detonation 
front"  and  "behind  the  detonation  front"  more  loosely,  in  terms  of  the 
measured  detonation  velocity  and  the  measured  velocities  of  shocks 
transmitted  through  normal  impact  on  various  inert  materials  in  con- 
tact with  the  explosive  (see,  for  example,  Deal4 5).  The  transmitted 
velocities  referred  to  are  those  that  would  be  found  in  a specified 
test  configuration,  by  extrapolating  back  to  the  explosive-inert 
material  interface,  with  the  effects  of  sharp  reaction  zone  spikes 
ignored.  Energy  releases  which  influence  one  or  more  of  these  measure- 
ments will  be  presumed  to  have  taken  place  at  the  detonation  front. 
Energy  releases  which  do  not  affect  the  detonation  velocity  or  any 
of  the  transmitted  shock  velocities  will  be  associated  with  flow  behind 
the  detonation  front.  In  the  same  vein,  it  will  be  assumed  that  the 
Chapman- Jouguet  (CJ)  condition  holds  at  the  detonation  pressure- 
particle  velocity  point  determined  by  the  transmitted  shock  experi- 
ments . 

A number  of  models  for  the  Amatex-20  detonation  have  been  sug- 
gested. Mader6 7  assumes  that  fifty  percent  of  the  AN  reacts  at  the 
detonation  front  and  that  the  remainder  remains  inert.  Thermo- 
chemical calculations,  assuming  complete  chemical  equilibrium  of  the 
products  from  the  elemental  composition  of  the  mixture  of  RDX , TNT, 
and  reacting  AN,  then  produce  the  observed  detonation  velocity  and  a 
reasonable  detonation  pressure.  Kamlet  has  proposed  a model  in 
which  the  AN  and  the  RDX/TNT  mixture  both  react  completely  but 
separately  at  the  detonation  front,  each  going  to  its  own  detonation 
products.  One  would  then  have  a heterogeneous  mixture  of  gas  pockets 
in  which  further  mechanical  and  chemical  interaction  takes  place  after 
the  detonation  front  has  passed.  Hershkowitz  and  Rigdon  point  out 
that  the  examination  of  wall  expansion  data  from  thin  and  thick  walled 
cylinders  shows  that  energy  added  to  the  system  due  to  the  AN  must 


4.  H.  M.  Sternberg,  "On  the  Mathematical  Theory  of  the  Chapman 
Jouguet  State",  Astronautica  Acta,  Vol.  15,  pp  359-369,  Pergamon 
Press,  New  York,  1970. 

5.  W.  E.  Deal,  "Measurement  of  the  Reflected  Shock  Hugoniot  and 
Isentrope  for  Explosion  Reaction  Products",  Phys . Fluids,  Vol.  I, 

No.  6,  Nov-Dec  1958.  Phys.  Fluids  1,  523  (1958). 

6.  C.  L.  Mader  A.  Popolato,  and  A.  W.  Campbell,  "EOS  of  Amatex" . 
Los  Alamos  Scientific  Laboratory,  Los  Alamos,  New  Mexico,  LA-5521-PR 
Joint  Services  Explosive  Program,  Quarterly  Report,  August  16 
through  November  30,  1973,  compiled  by  A.  Popolato. 

7.  Mortimer  J.  Kamlet.  Private  Communication. 
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appear  before  the  cylinder  walls  have  moved.  For  our  purpose  here 
this  means  that  the  additional  energy  must  be  added  at  the  detonation 
front  or  very  close  behind  it. 

The  computational  model  adopted  here  for  use  with  Amatex-20  is 
based  on  the  facts  that  its  density,  detonation  velocity,  and  deton- 
ation pressure  are  close  to  those  of  TNT,  but  the  wall  velocities 
found  in  cylinder  expansion  work  are  somewhat  greater  than  those 
obtained  with  TNT.  To  match  these  data  an  equation  of  state  of  the 
Jones-Wilkins-Lee  (JWL)  form8,  with  constants  chosen  to  give  the 
correct  detonation  velocity  and  an  estimated  detonation  pressure, 
was  combined  with  an  energy  addition  of  200  calories  per  gram  to  each 
computation  cell  one  microsecond  after  the  detonation  front  had 
passed  through  the  cell.  The  short,  one  microsecond,  delay  might 
correspond  to  an  induction  time.  This  model,  when  used  in  the 
explosion  dynamic  computations  to  be  described  here  gave  results  in 
good  agreement  with  the  existing  experimental  data.  However,  no 
claim  is  made  for  uniqueness.  Other  combinations  of  equations  of 
state,  equation  of  state  constants,  delay  times,  and  energy  releases, 
will  no  doubt  work  as  well. 

In  Section  II,  the  JWL  equation  of  state  constants  for  Amatex-20 
are  worked  out,  and  the  energy  addition  scheme  used  in  the  explosion 
dynamic  calculations  is  described.  Section  III  contains  the  existing 
cylinder  expansion  data,  the  calculated  cylinder  wall  velocities, 
and  a discussion  of  the  relatively  rapid  acceleration  of  the  cylinder 
wall  observed  at  late  times  in  the  LASL  work"  with  Amatex-20.  The 
two-space  dimension,  transient,  explosion  dynamic  computations  to 
get  the  cylinder  wall  velocities  were  made  with  a computer  code  which 
is  described  in  detail  in  a previous  report9.  This  is  a Lagrangian 
code  that  treats  the  cylinder  wall  as  a collection  of  mass  points, 
i.e.,  the  metal  mass  is  accounted  for  by  increasing  the  masses 
associated  with  the  outermost  grid  points.  Otherwise  the  logic  is 
similar  in  many  respects  to  that  of  the  HEMP  code10. 

The  results  of  the  transient  spherical  flow  calculations  for 
the  BLU61B  and  BLU63B  bomblets  are  in  Section  IV. 

II.  EQUATION  OF  STATE  SCHEME  FOR  AMATEX  20 

For  use  in  the  explosion  dynamics  computations  we  need  an  equation 
of  state  of  Amatex-20,  together  with  mechanisms  for  creating  the 

5"!  Brigitta  M.  Dobrat  z , Ed.,  Properties  of  Chemical  Explosives  and 
Explosive  Simulants,  Lawrence  Livermore  Laboratory,  University  of 
California/Livermore,  UCRL-51319,  15  December  1972. 

9.  H.  M.  Sternberg,  "Computation  of  Weight,  Velocity,  and  Angular 
Distributions  of  Fragments  from  Naturally  Fragmenting  Weapons" 

Naval  Ordnance  Laboratory,  White  Oak,  Silver  Spring,  Maryland, 

NOLTR  74-77,  17  July  1974. 

10.  Mark  L.  Wilkins,  "Calculation  of  Elastic-Plastic  Flow,"  in  methods 
of  Computational  Physics,  Vol.  3,  edited  by  B.  Alder,  S.  Fernbach, 
and  M.  Rotenberg,  Academic  Press,  New  York,  1964. 
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detonation  wave  and  for  releasing  additional  energy  due  to  chemical 
reaction  behind  the  front. 

The  JWL  form8  was  used  for  the  equation  of  state.  This  can  be 
written 


p = A ( 1 - R i V + B (1 

i 


co 


R V' 

2 


-R  V 

„ 2 


(1) 


where  p is  the  pressure,  E is  the 
V is  the  relative  specific  volume 
are  constants.  The  corresponding 
is 


internal  energy  per  original  cc, 

(=  v/v  ) , and  A,B,w,Ri,  and  R2, 
expansion  adiabat,  with  C constant, 


-R  V -R  V , . , . 

p = Ae  1 + Be  2 + CV  (a)  (2) 


The  constants  A,  B,  and  C give  p in  megabars. 

Reference  8 contains  a table  of  constants  for  various  explosives, 
for  use  in  the  JWL  equation.  For  most  common  military  explosives 
containing  only  carbon,  hydrogen,  oxygen,  and  nitrogen,  (CHNO)  R2  is 
between  0.9  and  1.2,  cu  is  about  0.35,  and  Rj  is  between  4.2  and 
4.6.  A choice  of  Rj  outside  of  a small  range,  which  is  different 
for  different  explosives,  will  cause  the  condition  that  the  sound 
speed  be  positive  to  be  violated  as  the  gas  expands.  In  general, 
increasing  the  value  of  Rj  inside  this  permissible  range,  with  a 
corresponding  change  in  A and  B to  keep  the  detonation  conditions 
fixed,  will  have  practically  no  effect  on  the  expansion  from  the  CJ 
state  to  V=l.  The  larger  value  of  Ri  will,  however,  cause  the  gas 
to  do  more  work  in  the  expansion  from  about  V=1  to  V=5,  and  less  for 
V > 5 . 

The  procedure  adopted  here  has  been  to  fix  Ri,  R2,  and  w,  and  to 
use  the  detonation  velocity,  the  detonation  pressure,  and  the  energy 
released  at  the  detonation  front  to  get  A,  B,  and  C in  Eq . (2) . This 

was  done  as  follows: 

Consider  a layer  of  gaseous  detonation  products  at  the  detonation 
front  and  its  behavior  when  the  time  is  too  short  for  later  energy 
additions  to  have  an  influence,  e.g. , when  a shock  is  transmitted 
into  a low  density  material.  Assume  that  the  extremely  rapid  expan- 
sion that  occurs  here  is  isentropic  and  that  it  starts  at  the  CJ 
state.  One  then  has,  neglecting  the  ambient  pressure  ahead  of  the 
detonation  front, 
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P.  = Po  Du, 

(3) 

D = u3  + c. 

(4) 

Vj  D = D - u,  = c. 

(5) 

(dp/ dV)^  = -po  c 2 /V/  = -po  D2 

(6) 

E,  = Q + P,  (1  - v,  )/2 

(7) 

Here  u is  the  particle  velocity,  c is  the  sound  speed, 
density  of  the  undetonated  explosive,  and  the  subscript 
the  CJ  condition.  In  Eq.  (7),  Q is  the  energy  released 
sive  at  the  detonation  front  (mbar-cc/original  cc) . It 
include  any  later  energy  addition. 

pQ  is  the 
J denotes 
by  the  explo 
does  not 

The  energy  Ej  in  Eq.  (6)  being  the  work  done  in  the 
expansion  of  the  gas  from  the  CJ  state  is 

00 

complete 

E =1  p dV 

(8) 

J 

v. 

Equation  (2),  together  with  its  derivative  and  integral,  and 
Eqs . (3) -(8) , now  lead  to  three  linear  equations  which  can  be  solved 

for  A,  B,  and  C,  once  pD  , D,  p.,  and  Q are  specified  (recall  that 
Rj,  R2  , and  go  are  chosen  in  advance).  These  equations  are: 


-R  V, 


A + e B + V 


-(«  + 1) 


C = P„ 


R e-RiVj  A + R e"R2V;r  B + (u>  + 1)  V ~(w  +2)  C = p D2  , 

12  jo 


(9) 

(10) 


-R  V 


-R  V. 


e i'J  , ^ e *Vr  - • V. 


—co 


A + 


B + 


C = Q + p,/  (2  PoD  ) 


(ID 


co 


The  density  pe  = 1.603  g/cc  and  the  detonation  velocity  D = 0.7030 
cm/iisec  were  used  for  Amatex  20  (see  Refs.  2 and  3)  in  Eqs.  (10)  and 
(11).  The  detonation  pressure  needed  in  Eqs.  (9)  and  (11)  was 
estimated  with  the  Kamlet  formula11  for  CHNO  explosives,  viz.. 


11 . Mortimer  J.  Kamlet  and  S.J.  Jacobs,  "Chemistry  of  Detonations, 

I.  A Simple  Method  for  Calculating  Detonation  Properties  of  C-H-N-0 
Explosives,"  J.  Chem.  Phys.  vol.  48,  pp  23-35,  1968. 
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P,  =1.53  [D/(1.30  + l/p„  ) ] 2 


(12) 


with  D in  cm/ysec,  pe  in  g/cc,  and  p^in  mbar. 

With  the  D and  p«,  values  given  above  the  value  p,  = 0.204  megabars 
was  calculated  from  Eq . (12).  The  constants  Ri , R2 , and  w were  taken 

as  4.6,  1.0,  and  0.35  respectively.  The  energy  release  at  the  deton- 
ation front  was  arbitrarily  set  at  Q = 0.0536  mbar-cc/original  cc 
(800  calories  per  gram).  The  solution  of  Eqs.  (9)-  (11)  was  then 

A = 4.9521  B = 0.08058  C = 0.001954. 

The  equation  of  state  Eq . (1)  was  augmented  with  a simple  scheme 

for  bringing  the  effects  of  the  chemical  reactions  into  the  explosion 
dynamic  computations.  Note  below  that  all  of  the  energy  that  will 
eventually  be  used  is  put  into  each  computation  cell  originally,  but 
only  part  of  this  energy  may  be  used  to  calculate  the  pressure, 
depending  upon  the  values  of  the  progress  variables  F]  and  F2  , 
which  are  initially  zero. 

If  Eq . (1)  is  written 


p = f (E,  V) 


(13) 


then  the  form  used  in  the  computations  is 


p = f (E,  V)  F 

1 K 


(14) 


where  F i*  is  a burn  fraction  corresponding  to  the  fraction  of  the 
computation  cell  considered  detonated  (see  the  description  of  the 
HEMP  code  1 0 ) , and 


E = E - (1  - F ) Q . 

2 k 2 


(15) 


F 


l 


F 

l 


The  variable  FiKis  calculated  from 
= max  |(D-T  - Z K>0  )/(ZK+1>o  - Z,  c ),  (1 

= min  | 1 , F /r  } . 


VK.„  >/<!  - V,  >,| 

(16) 
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where  Z is  the  axial  distance,  k is  the  Lagrangian  label  for  the 
grid  points  in  the  axial  direction,  and  T is  the  time  from  the  start 
of  the  detonation  at  Z = 0.  The  value  1.0  was  used  for  r. 

The  constant  Q2,  in  Eq . (15),  is  the  additional  energy  per  ori- 
ginal cc,  added  behind  the  detonation  front.  The  progress  variable 
F2  would,  in  more  complicated  cases,  be  found  by  integrating  the 
reaction  rate  dF2/dt,  which  would  be  a function  of  the  equation  of 
state  variables.  In  the  present  calculations  F2  remains  equal  to 
zero  until  a specified  number  of  microseconds  after  the  detonation 
front  has  passed  the  k ' th  cell.  At  this  time,  all  of  the  additional 
energy  02  is  made  available  for  the  pressure  calculation  by  abruptly 
changing  F2  from  zero  to  one. 

Note  above  that  the  total  energy  put  into  each  cell  at  the  outset 
is  Q + Q2.  However,  this  is  not  all  available  for  the  pressure 
calculation  until  both  Fj  and  F2  become  equal  to  one. 

In  the  hydrodynamic  computations,  Eqns.  (14)  and  (15)  are 
used  together  with 

dE  + (p  + q)  dV  = 0 

to  find  p and  E,  once  V is  known.  The  variable  q is  the  artificial 
viscosity . 

Table  I is  a summary  of  the  JWL  form  equation  of  state  and  the 
related  computational  scheme  worked' out  for  Amatex-20.  The 
additional  energy  Q2  was  taken  equal  to  200  calories  per  gram 
because  this  resulted  in  calculated  cylinder  wall  velocities 
which  matched  the  existing  data  (see  Section  III). 
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TABLE  I 

Equation  of  State  Scheme  Used  for  Amatex-20 
JWL  Equation  of  State 


P = A ( 1 

id  = 0.35 
Ri  = 4.6 
R2  = 1.0 


id  . - R V , _ , , 

) e i + B (1 


R V 

l 


A = 4 . 9521 
B = 0.08058 


ID 


R V' 


2 


R v , 
2 + 


ojE 


p0  = 1.603  g/cc,  D = 0.7030  cm/ysec,  VT  = 0.7425 
Eq  = Q + Q2  = 0.0675  mbar-cc/or iginal  cc  (1000  cal/g) 
Q = 0.0536  mbar-cc/original  cc  (800  cal/g) 

Q2  = 0.0139  mbar-cc/original  cc  (200  cal/g) 

For  detonation  wave  computation 


p = f (E,  V)  F 

1 K 


E = E - (1  - P ) Q 

2 k 2 


t = max  (D*T  - Z )/ (Z 

1 v*.  *<  + 1,  o 


(1  - v»,„  i/d-v,  )} 


F = min  {l,  E /r]  , r = 1.0 

IK  1 K 


0,  T < T K + 1 ysec 


F = 1 

2k 

'l,  T TK  = 1 ysec 


TK = time  detonation  front  reaches  the  end  of  the  k ' th  computation 
cell . 
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III.  CYLINDER  TESTS 

A great  deal  of  testing  has  been  done  with  a Lawrence  Livermore 
Laboratory  (LLL)  configuration  described  in  the  Fourth  Detonation 
Symposium  volume12.  This  consists  of  a one  inch  inside  diameter 
copper  tube,  12  inches  long,  with  a 0.1022  inch  wall  thickness.  The 
explosive  is  plane  detonated  at  one  end  and  the  position  of  the  out- 
side wall  versus  time  in  a plane  perpendicular  to  the  axis  and 
seven  inches  from  this  end  is  found  from  streak  camera  records.  This 
is  differentiated  to  get  the  outside  wall  velocity  in  the  observation 
plane . 

The  LLL  cylinder  test,  scaled  up  to  a four  inch  inside  diameter, 
was  used  at  both  PA2  and  LASL3  to  determine  the  metal  accelerating 
abilities  of  Amatex-20  and  Amatex-40  (40/40/20  RDX/TNT/AN)  relative 
to  a standard.  The  velocity  data  from  the  Amatex-20  loaded  cylin- 
ders, with  explosive/metal  mass  close  to  0.4,  will  be  used  here  to 
get  the  energy  to  be  added  behind  the  detonation  front  in  the  explo- 
sion dynamic  computations.  Before  doing  this  the  data  from  the 
different  laboratories  will  be  compared. 

The  variations  in  the  compositions  of  the  standards  used  at  the 
different  laboratories  are  insignificant.  The  standard  used  at  PA 
was  60/40  Cyclotol  (60/40  RDX/TNT) . The  60/40  refers  to  the  starting 
percentages,  by  weight,  in  the  formulation  process.  The  finished 
explosive  is  in  the  61/39  to  64/36  range.  The  standard  used  at  LASL 
and  LLL  was  Composition  B,  Grade  A.  Here  the  starting  materials  are 
59.5/39.5/1.0  RDX/TNT/Wax . At  LASL,  the  finished  explosive  is  about 
63/37/0.4  RDX/TNT/Wax11.  In  the  LLL  table  of  equation  of  state 
constants8  64/36  RDX/TNT  is  listed  for  Composition  B,  Grade  A. 

The  measured  solid  explosive  densities  were  1.694  g/cc  in  the  PA 
work2,  1.701  at  LASL3,  and  1.717  at  LLL12. 

Figure  1 contains  plots  of  the  outer  wall  velocities  reported 
by  the  three  laboratories  for  the  60/40  Cyclotol  and  Composition  B, 
Grade  A,  all  scaled  to  the  four  inch  inside  diameter  cylinder. 

Both  the  PA  and  the  LASL  measurements  were  made  in  a plane  24  inches 
from  the  booster  end  of  the  cylinder.  The  corresponding  scaled 
distance  for  the  LLL  work  is  28  inches.  Calculations  show  that  R-R 
versus  t at  the  two  positions  should  be  essentially  the  same.  Only 
two  of  the  LLL  velocities  are  shown  in  Fig.  1.  These  are  the 
velocities  listed  in  Ref.  12.  Note  in  Fig.  1,  that  at  R-Rq  = 76  mm 
(the  19  mm  in  the  LLL  test  scaled  to  a four  inch  I.D.  cylinder),  the 
LLL  and  LASL  velocities  are  practically  identical,  and  the  PA  velocity 
is  about  one  percent  greater.  At  R-R  = 104  mm  the  PA  velocity  is 
about  two  percent  greater  than  the  LA§L  velocity. 

12.  J.W.  Kury,  H.C.  Hornig,  E.L.  Lee,  J.L.  McDonnel,  D.L.  Ornellas, 

M.  Finger,  F.M.  Strange,  and  M.L.  Wilkins,  "Metal  Acceleration  of 
Chemical  Explosives,"  in  Proc . Fourth  Symp.  (International)  on 
Detonation,  Office  of  Naval  Research,  Rpt.  ACR-126,  U.S.  Gov't 
Printing  Office,  Washington,  D.C.  1965. 

13.  A.  Popolato.  Private  communication. 
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Outer  cylinder  wall  velocities  obtained  by  PA2  and  LASL 3 with 
Amatex-20,  in  the  LLL  test  scaled  to  the  four  inch  I.D.,  are  also 
shown  in  Fig.  1.  At  R-Rq  = 76  mm,  the  LASL  velocity  is  about  three 
percent  lower.  However,  between  R-R  = 90  mm  and  R-R  = 115  mm  the 
LASL  data  contains  a rapid  increase  in  velocity  which  does  not  appear 
in  the  PA  work. 

In  connection  with  the  possibility  of  very  late  energy  addition 
due  to  chemical  reaction,  it  is  of  interest  to  estimate  how  much 
energy  would  have  to  be  added  to  the  system  in  order  to  produce 
the  rapid  velocity  increase  observed  in  the  LASL  work. 

To  make  this  estimate^  assume  radial  flow  with  uniform  pressure 
p,  density  p,  and  internal  energy  E,  which  are  functions  of  time 
only.  Treat  the  metal  as  a tube  with  the  correct  mass,  but  no 
thickness.  Assume  also  that  the  gas  obeys  a V law,  and  expands 
adiabatically , in  the  time  interval  ti«t<t2,  where  the  sub- 
scripts 1 and  2 mark  the  ends  of  the  rapid  acceleration  period.  Let 
R be  the  radius  of  the  metal  cylinder  and  u the  velocity  ( = dR/dt) . 
The  equation  of  motion  can  then  be  written 


R°  l°  du  _ 

2C/M  dt  RP ' 


(17) 


where  pa  is  the  solid  explosive  density,  C/M  is  the  solid 
explosive  mass/metal  mass,  and  R is  the  initial  radius.  The 
gas  behavior  assumptions  imply 

P/P  = (P/P  )r  (13) 

1 1 

and  the  conservation  of  mass  requires 

p/p  = R2  / R2 . (19) 

l l 


From  (17)  - (19), 


i - 2 r 

I<  dZ/Dt  = p r , 

l 

where  K = R0  pc  / ( 2C/M)  , r = R/R  and  Z = dr/dt. 


(20) 
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Equation  (20)  has  the  integral 


K 


1 


1 


(i-r)  j 


(21) 


Suppose  that  additional  internal  energy  is  added  abruptly  at  time  t . 

This  will  increase  the  pressure  but  leave  the  velocity  and  density  1 

at  t=t  unchanged.  Call  the  quantities  changed  by  the  energy  addition 

E , p j u , r . The  y law  assumption  means 
112  2 


E = 

1 p~cr-T) 

i 


a nd  E = 

1 p (r-1) 

i 


From  Equation  (21), 


R 


[u2  - u2]  = [1  - r 2 (l"  Y)  ] . 

~2  2 1 t ~ L 


1 


At  the  same  radial  expansion  (r2=r2),  using  (21)- (23), 

/\ 

E 


u2  - u2 


1 _ 2 


u2  - u2 
2 1 


(22) 


(23) 


(24) 


In  Fig.  1 the  LASL  wall  velocity  of  the  Amatex-20  loaded  cylinder 
is  1.485  mm/psec  (this  is  Ui)  at  R-Ro=90  mm.  The  velocity  U2  at 
R-R0=115 . 6 mm  is  1.560  mm/psec.  The  extrapolated  velocity  at 
R-R0=115.6  mm,  taken  from  the  dashed  section  in  Fig.  1,  is  1.510 
mm/psec.  This  is  U2.  With  the  above  values  for  ui,  U2,  and 
u 2 , we  get,  with  Eq . (24), 


E /E  = 3.05.  (25) 

1 1 

Transient  explosion  dynamic  computations  in  two- space  dimensions 
show  that  at  R-Ro=90  mm  about  15%  of  the  chemical  energy  released  in 
the  early  stages  remains  as  internal  energy  of  the  gas.  Using  this 
value,  together  with  the  energy  release  in  Table  I (1000  cal/g)  and 
Eq . (25),  shows  that  about  300  additional  calories  per  gram  would  have 

to  be  added  at  R-Ro=90  mm  in  order  to  produce  the  velocity  increase 
observed  in  the  R-Ro=90  to  115  mm  range. 

The  computer  program,  MASSPT  described  in  Ref.  (9)  was  used  for 
the  explosion  dynamic  computations  involving  the  cylinder  test 
configuration.  This  program  treats  the  metal  as  a tube  with  mass 
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but  no  thickness,  i.e.,  as  a set  of  grid  points  with  extra  mass, 
whose  motions  are  found  along  with  the  gas  flow.  The  gas  dynamics 
following  detonation  is  Lagrangian,  with  provision  for  slippage 
along  the  metal  boundary.  Additional  routines,  not  needed  in  the 
weapon  calculations  described  in  Ref.  (9),  were  added  in  order  to 
permit  the  gas  to  flow  out  of  the  ends  of  the  cylinder.  In  the 
computation  the  cylinders  were  4 inches  in  inside  diameter,  38  inches 
long  (36  inches  plus  2 inches  for  the  plane  wave  initiator) , and 
plane  initiated  at  one  end.  A 100X6  computation  grid  was  used, 
making  the  computation  cells  0.847  cm  in  the  radial  direction  and 
0.965  cm  in  the  axial  direction.  To  compensate  for  the  coarse 
zoning  the  energy  was  released  more  quickly  than  in  the  HEMP 
code10,  by  taking  r in  Eq . (16)  equal  to  1.0,  instead  of  the 

2.5  value  suggested  in  Ref.  (10). 

The  MASSPT  code  gives  the  motion  of  a mass  tube  of  zero  thickness. 
The  motion  of  the  outside  wall  must  be  gotten  from  this  for  comparison 
with  the  experimental  data.  Let  R be  the  outside  radius,  R the 
inside  radius,  and  use  the  zero  subscript  for  the  initial  conditions. 
Assume  the  metal  is  incompressible.  Then  the  outside  and  inside  wall 
radii  and  velocities  are  related  by 


R2 


R' 


r; 


and,  using  Eq . (26), 


dR  _ R dR 

dt  ~ r ar ■ 


(26) 


(27) 


Also,  if  r is  a variable  radius  between  the  inner  and  out  walls, 


R2 


Rj 


and 


dR 

dt 


r dr 
R dt* 


(28) 


The  outer  wall  velocity  in  the  observation  plane  was  found  here  from 
the  calculated  mass  point  velocity  by  assuming  that  the  kinetic 
energy  of  a section  of  the  tube  is  equal  to  the  kinetic  energy  of  the 
corresponding  mass  point,  counting  only  the  radial  velocity.  We  then 
have,  with  v the  calculated  radial  velocity  of  the  mass  point,  and 
the  metal  assumed  incompressible, 


(29) 


I 
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Equations  (28)  and  (29)  lead  to  the  outer  wall  velocity 

r ~ , -i1'2 

dR  |l  - (R/R) 

dt  J?,n  ( R/£)  ' 


(30) 


where  R is  the  radial  position  of  the  mass  point  and  R and  R are 
related  by  Eq . (26).  A very  close  approximation  to  Eq . (30), 
in  the  range  of  interest,  is 


dR  _ 1 
dt  " 7 


[■ 


(31) 


The  curves  in  Fig.  1 are  plots  of  velocity  versus  distance 
calculated  with  the  MASSPT  program  for  the  plane  perpendicular  to  the 
axis  of  the  copper  cylinder  and  26  inches  from  the  initiation  end. 
Twenty  six  inches,  instead  of  24,  were  used  here  in  order  to  take 
the  plane  wave  booster  into  account.  The  calculated  mass  point  velo- 
cities were  corrected  with  Eq . (30)  to  get  the  plotted  outer  wall 

velocities.  The  use  of  this  correction,  which  implies  a uniform 
metal  density  and  instantaneous  communication  between  the  inner  and 
outer  cylinder  walls,  probably  accounts  for  the  calculated  outer 
wall  velocities  being  a bit  higher  than  the  experimental  values 
during  the  early  motion. 

The  calculations  for  the  cylinders  containing  TNT  and  60/40 
Cyclotol  were  made  with  the  JWL  equation  of  state  and  the  JWL  equation 
of  state  constants  listed  in  Ref.  (8).  The  Composition  B,  Grade  A 
constants  in  Ref.  (8)  and  the  solid  explosive  density  p=1.701  g/cc 
was  used  to  model  both  the  PA  and  LASL  experiments.  This  density 
is  the  experimental  value  reported  for  the  LASL  charges  . 

The  calculated  results  shown  for  Amatex-20  in  Fig.  1 were  obtained 
with  the  equation  of  state  scheme  in  Table  I.  This  has  the  JWL  form 
with  the  constants  found  by  the  method  in  Sec.  II.  Also  (see  Table  I), 
in  order  to  match  the  experimental  data,  200  cal/g  were  added  to 
each  computation  cell  one  microsecond  after  the  detonation  front 
passed  through  it.  This  choice  produces  velocities  at  R-R0=76  mm  that 
are  roughly  midway  between  the  LASL  and  the  PA  values. 

If  one  wanted  a computation  scheme  that  would  produce  only  the 
LASL  results  for  Amatex  20  in  the  4 inch  I.D.  cylinder  including  the 
very  late  acceleration  (Fig.  1) , then  an  extension  of  Eq . (15)  could 

be  used,  namely, 

E = E - (l-F)Q  - (1  - F ) 0 • 

2 2 3 3 
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Here,  Q2=150cal/g,  03=300  cal/g,  F2  and  F3  are  zero  initially, 

F 2 is  set  equal  to  one  after  a one  microsecond  delay,  and  F3  is 
set  equal  to  one  when  R-Ro=90  mm  in  the  cross  section  intersecting 
the  computation  cell.  The  total  energy  release,  including  the 
800  cal/g  at  the  detonation  front,  would  then  be  1250  cal/g.  This 
is  reasonable  for  the  elemental  composition  of  the  solid  explosive. 
However,  a rationale  for  the  very  late  energy  release  is  missing. 

IV  COMPUTATIONS  FOR  BOMBLETS 

The  BLU-61B  and  BLU-63B  bomblets  are  small  spherical, steel  shells 
(consisting  of  two  hemispherical  shells  fastened  together)  containing 
a grenade  fuze,  a booster  and  a high  explosive  fill.  The  fragment 
size  is  determined  by  scoring  the  metal  surface.  The  fuze  initiates 
the  booster  at  an  off  axis  spot  about  one  half  inch  from  the  center. 

Transient,  one  space  dimension,  spherical  flow  calculations  were 
made  with  a Lagrangian,  artificial  viscosity,  explosion  dynamics 
code.  Fifty  computation  cells  were  used.  The  innermost  grid  point 
was  fixed  at  1.27  cm  from  the  center  for  the  entire  computation,  i.e., 
the  fuze  was  modeled  by  a rigid  sphere  with  a one  half  inch  radius, 
and  a spherical  detonation  wave  was  initiated  on  the  surface  of  this 
sphere.  The  artificial  viscosity  was  a sum  of  a linear  and  a qua- 
dratic function  of  the  particle  velocity.  The  metal  casing  was 
modeled  by  adding  its  mass  to  the  mass  associated  with  the  outer- 
most grid  point.  Dimensions  used  in  the  computations  for  the 
BLU-61B  and  the  BLU-63B  munitions  are  shown  in  Fig.  2.  The  equations 
of  state  used  for  the  cylinder  expansion  calculations  described  in 
Sec.  Ill,  were  used  here  namely,  the  scheme  for  Amatex-20  listed 
in  Table  I,  and  the  JWL  equation  for  Composition  B,  Grade  A listed 
in  Ref.  (8)  (with  pQ=l . 701g/cc  and  the  corresponding  Eo=0.0842  mbar- 
cc/or ig . cc ) . 

The  hemispheres  of  the  submunitions  move  apart  initially  before 
they  break  up.  This  results  in  some  venting  of  the  product  gas 
and  a rarefaction  wave  which  moves  inward  and  effects  the  metal 
velocity  near  the  edges  of  the  hemispheres.  This  rarefaction  and 
the  related  pressure  drop  are  ignored  in  the  transient  one  soace 
dimension  calculations. 

Figures  3 and  4 contain  plots  of  the  calculated  casing  velocity 
versus  inside  radius.  Note  that  the  velocity  keeps  increasing 
because  there  is  nothing  in  the  computation  to  account  for  the  fragmen- 
tation of  the  casing.  To  get  a calculated  fragment  velocity,  one 
must  decide,  usually  from  photographic  evidence,  upon  the  radius  at 
which  the  casing  breaks  up  into  fragments  and  acceleration  essen- 
tially stops.  This  radius  is  not  available  here.  We  can,  however, 
try  to  work  back  from  the  observed  fragment  velocities  in  order  to 
estimate  the  radius  where  the  casing  breaks  up. 
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The  initiation  spot,  in  the  grenade  fuze,  is  approximately 
one  half  inch  from  the  center  and  70  degrees  from  the  pole  of  one  of 
the  hemispheres.  Fragment  velocities  were  estimated  by  Ramsay14  at 
LASL,  from  flash  radiographs.  For  Amatex  20  loaded  bomblets  the 
estimated  velocities  at  the  70  degree  position  were  1.37±0.05  mm/ysec 
(4495±164  ft/sec)  for  the  BLU-61B  munition  and  1.1±0.05  mm/ysec 
(3609±164  ft/sec) for  the  BLU-63B  munition.  Figs.  3 and  4 show 
that  these  velocities  correspond  to  acceleration  cessation,  due  to 
casing  break  up,  at  46%  radial  expansion  (5.2  cm  inner  radius)  for 
the  BLU-61B  munition  and  at  23%  radial  expansion  (3.2  cm  inner  radius) 
for  the  BLU-63B  munition.  The  0.05  mm/ysec  standard  deviation 
in  the  measured  velocity  corresponds  to  a 3-4%  change  in  the  radial 
expansion  when  the  casing  breaks  up. 

In  the  cylinder  test  experiments2*3  it  was  shown  that  the  substi- 
tution of  40  parts  of  AN  for  RDX  in  60/40  Cyclotol , to  make  Amatex -20 , 
resulted  in  a kinetic  energy  imparted  to  the  metal  ratio  of  0.83-0.85. 
This  ratio  can  be  estimated  for  the  munitions  considered  here,  from 
Figs.  3 and  4,  by  assuming  that  the  kinetic  energy  ratio  is  the  square 
of  the  calculated  velocity  ratio.  This  leads  to  kinetic  energy 
ratios  of  0.86-0.88  for  the  BLU-61B  munition  and  0.81-0.87  for  the 
BLU-63B  munition.  The  exact  values,  in  these  ranges,  depend  upon  the 
radial  expansion  when  acceleration  ceases.  If  the  dashed  vertical 
lines  in  Figs.  3 and  4 are  used  (46%  radial  expansion  for  the 
BLU-61B  and  23%  for  the  BLU-63B) , the  calculated  ratios  are  0.88 
for  the  BLU-61B  and  0.81  for  the  BLU-63B  munitions. 

V REMARK  ON  FRAGMENTATION 

Amatex-20  has  been  found  unsuitable  for  use  in  the  BLU-61B  and 
the  BLU-63B  munitions  because  the  casings  do  not  break  along  some 
of  the  scored  lines  and  significantly  more  multiple  fragments  result 
than  are  obtained  with  the  preferred  fill  (60/40  Cyclotol) . At  one 
time  it  was  thought  that  this  could  be  explained  in  terms  of  the 
small  size  of  the  systems  and  the  run  distance  needed  to  get  the 
detonation  fully  started.  This  does  not  appear  to  be  the  case 
(see  Part  B of  this  report). 

A plausible  explanation  for  the  greater  number  of  multiple  frag- 
ments found  when  Amatex-20  is  used  lies  in  the  fact  that  the  ability 
of  this  explosive  to  fragment  metal  is  inferior  to  that  of  the  60/40 
Cyclotol.  The  density,  detonation  velocity,  and  detonation  pressure 
of  Amatex-20,  which  affect  the  fragmentation  are  very  much  like  those 
of  TNT.  The  values  of  these  quantities  that  we  have  been  using  here 
for  Amatex-20,  TNT,  and  60/40  Cyclotol  (see  Sec.  II  and  Ref.  8)  are 
listed  in  Table  II. 


14.  J.B.  Ramsay,  Bomb  Performance  Tests,  in  "Feasibility  of  Using 
Amatex  in  Small  Bombs,  Progress  Report  Sep  1 through  Oct  31,  1974, 
AFATL  project  order  4-065Q-11"  University  of  California,  Los  Alamos 
Scientific  Laboratory,  WX-3,  Dec  6,  1974,  letter  report  compiled 
by  A.  Popolato. 
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Table  II 

Densities,  Detonation  Velocities,  and  Detonation  Pressures 


of 

Amatex  20,  TNT, 

and 

60/40  Cyclotol 

Amatex  20 

TNT 

60/40  Cyclotol 

p0  (g/cc) 

1 . 603 

1 . 630 

1 .701 

D (cra/ysec 

) 0.703 

0.693 

0.798 

p (mbar) 

0 .204 

0.210 

0 .295 

Reference  (15)  lists  the  relative  fragmentation  abilities  of  various 
explosives  in  terms  of  the  average  weight  of  fragments  M,  weighing  more 
than  one  grain,  recovered  in  a particular  test  configuration.  Values 
of  M obtained  there  are  17  grains  for  cast  TNT  and  12  grains  for 
Composition  B,  which,  except  for  a fraction  of  one  percent  of  wax,  is 
identical  to  60/40  Cyclotol.  From  Table  II  we  can  expect  that  the 
Amatex-20  will  have  about  the  same  ability  to  fragment  metal  as  TNT. 

We  can  therefore  assign  values  of  M equal  to  17  grains  for  Amatex-20 
and  12  grains  for  the  preferred  fill  (60/40  RDX/TNT) . This  indicates 
that  the  metal  fragmenting  ability  of  the  preferred  fill  is  substan- 
tially better  than  that  of  the  Amatex-20.  It  is  most  likely  the 
reason  for  the  greater  number  of  multiple  fragments  observed  with 
Amatex-20  than  with  60/40  Cyclotol. 

VI  CONCLUSIONS 

An  equation  of  state  scheme  for  Amatex-20  which  combines  a JWL 
equation  with  an  energy  addition  of  200  cal/g  one  microsecond  after 
the  detonation  front  has  passed  can  be  used  in  hydrodynamic  calculations 
with  satisfactory  results. 

Cylinder  test  wall  velocities  from  LASL,  PA,  and  LLL  agree  within 

2 percent  for  60/40  Cyclotol,  for  expansions  to  twice  the  initial 

diameter.  For  Amatex-20  in  four  inch  I.D.  cylinders  the  measured  PA 

wall  velocity  is  about  3 percent  greater  than  the  LASL  velocity  at 

R-R  =76  mm. 
o 

A rough  estimate  shows  that  an  abrupt  addition  of  about  300  cal/g 

to  the  detonation  products  of  Amatex-20  is  needed  to  produce  a wall 

velocity  increase  observed  between  R-R  = 90  mm  and  R-R  = 115  mm  in 
J o o 

the  LASL  four  inch  I.D.  cylinder  tests. 


15.  H.M.  Sternberg,  "Fragment  Weight  Distributions  from  Naturally 
Fragmenting  Cylinders  Loaded  with  Various  Explosives"  Naval 
Ordnance  Laboratory,  White  Oak.  Silver  Spring,  Maryland  20910, 
NOLTR  73-83,  12  October  1973. 
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Results  of  hydrodynamic  computations  combined  with  bomblet 
fragment  velocities  measured  at  LASL  from  flash  radiographs  indicate 
that  casing  rupture  to  form  fragments  occurs  after  46  percent  radial 
expansion  in  the  BLU-61B  and  after  23  percent  radial  expansion  in  the 
BLU-63B  bomblets . 

The  failure  of  Amatex-20  loaded  bomblets  to  break  on  some  of  the 
scored  lines  is  probably  due  to  the  inferior  metal  fragmenting  ability 
of  this  explosive,  relative  to  the  preferred  fill. 


20 


NSWC/WOL/TR  75-94 
PART  B 

EXPERIMENTAL  STUDY  OF  AMATEX-20  INITIATION 
AND  DETONATION 

N.  L.  COLEBURN  and  H.  D.  JONES 


I . INTRODUCTION 

Measurements  of  the  initiation  and  detonation  of  Amatex-20 
(20/40/40  RDX,  TNT,  AN),  using  standard  boostering  systems, 

were  made  to  give  some  insight  into  the  fragmentation  results 
obtained  when  this  explosive  is  used  as  a fill  in  the  BLU-61B 
and  the  BLU-63B  munitions  (bomblets) . In  fragmentation  tests 
the  munition  casings  failed  to  break  along  some  of  the  scored 
lines,  with  resulting  multiple  fragments.  It  was  felt  that 
this  might  be  connected  with  the  response  of  the  explosive 
to  the  type  and  strength  of  the  boostering,  detonation  propa- 
gation from  asymmetrical  initiation,  or  AN  particle  size 
effects.  The  relative  insensitivity  and  lower  energy  output 
of  Amatex- 20 , compared  with  the  preferred  fill,  were  also 
possible  causes  for  the  inadequate  fragmentation  performance. 

II.  EFFECT  OF  DETONATOR  BASE  CHARGE  WEIGHT  ON  THE 

INITIATION  OF  AMATEX-20  BY  STANDARD  BOOSTERS 

In  the  first  phase  of  the  initiation  study,  the  aquarium 
technique  and  high  speed  streak  camera  photography  were  used 
to  study  the  effect  of  variations  in  the  strength  of  the 
shock  from  a detonator.  Mark  70  detonators,  containing  a 
base  charge  of  80  mg  of  PETN , and  Mk  71  detonators,  contain- 
ing a base  charge  of  40  mg  of  PETN,  were  used.  The  tests 
consisted  of  measuring  the  lateral  and  axial  shock  output 
from  A-5  boosters  in  an  experimental  setup  similar  to  the 
boostering  arrangement  in  the  bomblets.  The  boosters  were 
cylindrical  pellets  with  an  L/D  of  0.3.  Two  pellet  sizes 
were  used,  15.3  mm  in  diameter  by  5.2  mm  long,  and  13.0  mm 
in  diameter  by  4.2  mm  long.  The  average  densities  were 
1.70  g/cc  and  1.62  g/cc,  respectively. 

The  charge  arrangement  is  shown  in  Figure  5A.  The 
detonator  was  inserted  into  an  aluminum  holder  and  was  held 
in  contact  with  one  face  of  the  boostering  pellet.  The 
detonator  was  offset  from  the  center  of  the  booster  face 
so  that  intiation  occurred  first  at  a point  3.1  mm  from 
one  edge.  The  arrangement  was  positioned  in  a Plexiglas - 
walled  aquarium  and  aligned  with  the  smear  camera  slit 
at  90  degrees  with  respect  to  the  booster  side  and  at  the 
midpoint  of  the  booster. 
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Figure  5B  is  a typical  record  from  the  experiment 
where  a 12.7  mm  diameter  A-5  booster  was  used.  The  camera 
writing  speed  was  'Mmm/visec.  The  first  event  seen  in 
Figure  5B  is  the  detonation  wave  arrival  on  the  side  of 
the  booster  as  viewed  through  the  camera  slit.  The  record 
shows,  as  expected,  earliest  time  arrival  for  the  wave 
emerging  from  the  booster  edge  which  was  closest  to  the 
point  of  initiation.  The  maximum  arrival  time  difference 
between  detonation  breakout  on  the  sides  of  the  cylindrical 
booster  was  vO . 4 microsec.  The  lateral  shock  traces  in 
the  water  also  are  clearly  delineated  in  Figure  5B.  Note 
the  large  differences  in  the  initial  decay  of  the  shock 
wave  propagation  from  each  side  of  the  booster. 

Shock-wave  velocities  were  obtained  by  differentiating 
the  distance- time  data  derived  from  microcomparator  readings 
of  the  shock  traces  in  Figure  5B.  The  velocities  were  con- 
verted to  peak  pressures  from  knowledge  of  the  shock  Hugoniot 
for  water1  . Figures  6 and  7 show  the  lateral  peak  pressure 
outputs  and  shock  wave  decay  derived  from  the  measurements . 

The  curves  in  Figure  6 were  obtained  from  measurements  of 
the  weaker  lateral  shock  (upper  trace  in  Figure  5B)  which 
propagates  from  the  side  of  the  cylindrical  booster  closest 
to  the  point  of  initiation.  Similarly  the  "strong  shock" 
curves  of  Figure  7 are  the  results  of  measurements  of  the 
shock  wave  which  originated  when  the  detonation  wave  reached 
the  side  farthest  from  the  point  of  initiation.  The  peak 
pressure  results  plotted  in  Figure  2 show  little  difference 
in  output  from  the  Mk  70  and  Mk  71  detonators.  In  the 
"strong  shock"  graph.  Figure  7,  an  anomaly  results.  The 
lateral  output  from  the  Mk  71  detonator  is  greater  than 
the  output  from  the  Mk  70  detonator  which  contains  the 
larger  base  charge  weight  of  PETN.  This  result  may  be 
due  to  the  experimental  error  in  the  measurements  and  the 
difference  in  the  densities  of  the  booster  charges.  The 
density  of  the  charge  fired  with  the  Mk  71  detonator  was 
1.819  g/cm3 while  the  booster  initiated  by  the  Mk  70  detonator 
had  a density  of  1.758  g/cm3.  The  "strong"  shock  output 
results  from  the  detonation  propagating  through  a greater 
thickness  of  explosive  and  thereby  would  be  more  affected 
by  the  density  difference  than  the  "weak"  lateral  shock 
output. 

III.  DETONATION  WAVE  ARRIVAL  MEASUREMENTS 

The  smear  camera  arrangement  shown  in  Figure  8 was  used 
to  observe  the  detonation  wave  arrival  in  Amatex-20. 
Measurements  were  made  of  wave  arrivals  on  the  charge  sides 
and  on  the  charge  face  opposite  the  detonator.  As  shown  in 
the  figure,  an  A-5  booster  was  imbedded  flush  with  one  face 
of  the  cylindrical  Amatex-20  charge.  Charges  containing 
500  micron  AN  and  charges  containing  10  micron  AN  were 
fired.  Fine  scratch  marks  were  cut  ^1  mm  deep  on  the  charge 

(16)  M.  H.  Rice  and  J.  M.  Walsh,  J.  Chem.  Phys.  2_6  , 824  (1957) 

(17)  N.  L.  Coleburn,  "The  Chapman- Jouguet  of  Several  Pure 
and  Mixed  Explosives",  NOLTR  64-58,  25  June  1964 
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surface  at  1.27  cm  intervals,  and  covered  with  scotch 
tape.  Detonation  wave  arrival  at  these  points  produced 
the  sharp  pips  noted  in  the  smear  camera  trace  of  Figure  9. 

The  detonation  wave  arrival  profiles  on  the  faces  of  the 
charges  were  somewhat  different  for  charges  containing  the  two 
different  ammonium  nitrate  particle  sizes.  The  wave  front  from 
the  500  micron  AN  charges  first  arrived  on  the  face  at  the  center. 
The  wave  front  from  the  10  micron  AN  charges  arrived  on  the  face 
at  the  point  offset  about  3 mm  from  the  center.  In  both 
cases,  subsequent  arrivals  of  the  wave  front  were  symmetric 
about  a line  perpendicular  to  the  face  through  the  point 
of  first  arrival  (see  Figure  10) . Approximate  detonation 
velocities  determined  from  the  breakout  on  the  charge  sides 
(see  the  top  and  bottom  traces  in  Figure  9)  are  a few 
percent  greater  for  the  charges  containing  10  micron  AN. 

Also  measurements  of  the  smear  camera  traces  in  Figure  9 
showed  that  the  velocity  of  the  circular  wave  front  on  the 
charge  face  is  a bit  greater  for  the  10  micron  AN  charge. 


IV.  DETONATION  PARAMETER  MEASUREMENTS 

Measurements  of  the  detonation  velocity  and  detonation 
pressure  of  Amatex-20  were  made  using  cylindrical  charges 
containing  AN  with  500  micron  particle  size.  The  detonation 
velocity  was  measured  directly  and  the  detonation  pressure 
was  inferred  from  shock  wave  measurements  using  the  aquarium 
technique  described  in  reference  (17) . In  these  measurements 
all  charges  were  eight  inches  long  and  two  inches  in  diameter. 
Boostering  was  accomplished  by  a one-inch  high,  two-inch 
diameter  pentolite-baratol  plane  wave  booster.  Wave  curva- 
ture effects  on  the  detonation  measurements  due  to  an  inade- 
quate L/D  (length  to  diameter)  were  minimized  by  the  above 
initiation.  The  charges  were  confined  in  2.00-inch  diameter 
copper  cylinders  with  0.233-inch  thick  walls.  Holes  of 
0.1-inch  diameter  were  drilled  through  the  cylinder  wall  on 
a line  extending  down  the  cylinder  length  at  0.5-inch  intervals. 
The  holes  were  filled  with  pressed  RDX  and  covered  with 
scotch  tape.  By  aligning  the  smear  camera  slit  to  bisect 
the  holes  and  to  extend  beyond  the  end  of  the  charge,  the 
detonation  velocity  of  the  explosive  and  the  velocity  of 
the  shock  transmitted  into  the  water  from  the  end  of  the 
charge  were  measured.  Cast  charges  with  densities  of 
1.603  - 0.002  g/cm3  were  used.  The  detonation  velocity 
measured  was  6790  i 20  m/sec.  The  initial  shock  velocity 
in  the  water  was  5440  - 35  m/sec.  The  corresponding  pressure, 

PH  0 and  particle  velocity  uH  are  121  kbar  and  2250  m/sec 
2 2 
respectively.  The  detonation  pressure  PT,  is  given  by 

solving  by  iteration  the  equations 
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P = P 

J ph20 


1 - 


(k  - 1)uH20  - (k  -1)D 


2kD 


P = 


poD 

k+1 


2k 

k-1 


/ 


where  k is  the  isentropic  exponent,  D the  detonation  velocity 
and  p , the  charge  density.  Substituting  the  above  values 

in  these  equations  yields  Pj  = 179  kbar  and  k = 3.12  for 
cast  Amatex-20  with  500  micron  AN  particle  size.  The 
densities  of  the  charge  in  the  tests  are  about  1%  lower  than 
the  charges  used  at  LASL,  and  assuming  an  error  of  ^1%  in 
measuring  the  detonation  velocity  and  shock  velocity,  an 
analysis  of  equations  and  indicates  the  detonation  pressure 
of  179  kbar  may  be  too  low.  This  measured  value  of 

the  detonation  pressure  is  ^10%  lower  than  the  pressure 
of  201  kbar  which  was  estimated  at  LASL  from  work  with 
larger  charges. 

V.  CONCLUSIONS 

1.  Amatex-20,  with  either  500  micron  or  10  micron 
particle  size  AN, can  be  initiated  satisfactorily  with  a 
standard  A- 5 booster. 

2.  No  differences  in  the  initiation  of  the  A-5  booster 
are  found  when  detonator  base  charge  weights  of  40  and  80 
mg  of  PETN , respectively,  are  used. 

3.  The  measured  detonation  velocity  of  Amatex-20 
charges  containing  500  micron  AN  was  6790  - 20  m/sec. 
Detonation  velocities  inferred  from  wave  profile  measure- 
ments, in  the  wave  arrival  experiment,  were  slightly 
greater  for  the  charges  containing  10  micron  AN  than  for 
those  containing  500  micron  AN. 

4.  The  measured  detonation  pressure  of  Amatex-20  con- 
taining 500  micron  AN  was  179  kbar,  for  a cast  charge  with 

a density  of  1.60  g/cc.  This  value,  obtained  in  the  aquarium 
experiment,  is  about  10%  smaller  than  that  reported  by  LASL, 
where  larger  charges  were  used. 


(18)  University  of  California,  Los  Alamos  Scientific  Laboratory, 
Los  Alamos,  New  Mexico,  "Feasibility  of  Using  Amatex  in  Small 
Bombs,  Progress  Report  June  16  through  August  30,  1974",  AFTL 
Project  Order  4,  065Q-111. 
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FIG.  1 CALCULATED  AND  EXPERIMENTAL  OUTER  WALL  VELOCITIES  OBTAINED  WITH  LLL  CYLINDER  TEST  SCALED  UP 
TO  4 INCH  INSIDE  DIAMETER.  RESULTS  FOR  60/40  CYCLOTOL,  AMATEX-20,  AND  TNT  ARE  COMPARED. 
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BLU-61B 


-1.27  cm 


H 


-3.556  cm- 


BLU-63B 


FIXED  WALL 


METAL  CASING  PLUS  COVER 
MODELED  BY  ADDING  312.5 
GRAMS  TO  MASS  OF  OUTERMOST 
GRID  POINT 


DETONATION  PRODUCT  GAS- 
50  COMPUTATION  CELLS 


-1 .27  cm- 


■2.593  cm- 


FIG.  2 DIMENSIONS  USED  FOR  SPHERICAL  FLOW  COMPUTATIONS. 
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FIG.  3 BLU-61B  SUBMUNITION.  CALCULATED  CASING  VELOCITY  AT  70°  POSITION  VS 
INSIDE  RADIUS,  FOR  60/40  CYCLOTOL  AND  AMAT EX-20  FILLS. 
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INNER  RADIUS  (cm) 

FIG.  4 BLU  - 63B  SUBMUNITION.  CALCULATED  CASING  VELOCITY  AT  70° 

POSITION  VS  INSIDE  RADIUS,  FOR  60/40  CYCLOTOL  AND  AMATEX-20  FILLS. 
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FIG.  5A 


ARRANGEMENT  FOR  LATERAL  SHOCK  OUTPUT  MEASUREMENTS 
FROM  A-5  BOOSTERS  VS  DETONATOR  BASE  CHARGE  WEIGHT. 
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FIG.  5B  SMEAR  CAMERA  RECORD  OF  LATERAL  SHOCK  WAVES  IN  WATER 


FROM  DETONATION  OF  AN  A-5  BOOSTER. 
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FIG.  6 LATERAL  PEAK  PRESSURE  VS  DISTANCE  IN  WATER  FROM 
THE  PERIPHERY  OF  A-5  BOOSTERS.  (WEAK  SHOCK) 
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PEAK  PRESSURE  (KBAR) 
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FIG.  7 LATERAL  PEAK  PRESSURE  VS  DISTANCE  IN  WATER  FROM 
THE  PERIPHERY  OF  A- 5 BOOSTERS.  (STRONG  SHOCK) 
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SMEAR  CAMERA 
SLIT  IMAGE 


I 


B.  CAMERA  VIEW 


FIG.  8 ARRANGEMENT  FOR  DETONATION  WAVE  ARRIVAL  MEASUREMENTS 
ON  THE  PERIPHERY  AND  FACE  OF  AMATEX-20  CYLINDER  CHARGE 
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FIG.  9 SMEAR  CAMERA  RECORD  OF  DETONATION  WAVE  ARRIVALON  THE 
PERIPHERY  AND  FACE  OF  A CYLINDRICAL  AMATEX  20  CHARGE 
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FIG.  10  DETONATION  WAVE  ARRIVAL  TIME  DIFFERENCES  AT  POINTS  ON  THE  FACE  OF  AN 
AMATEX-20  CYLINDER  INITIATION:  MK  71  DETONATOR,  A-5  BOOSTER. 
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